Brightening of an Accretion Disk Due to Viscous Dissipation of 
Gravitational Waves During the Coalescence of Supermassive Black Holes 
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Mergers of supermassive black hole binaries release peak power of up to ~ 10 57 erg s _1 in gravi- 
tational waves (GWs). As the GWs propagate through ambient gas, they induce shear and a small 
fraction of their power is dissipated through viscosity. The dissipated heat appears as electromag- 
netic (EM) radiation, providing a prompt EM counterpart to the GW signal. For thin accretion 
disks, the GW heating rate exceeds the accretion power at distances farther than ~ 10 3 Schwarzschild 
radii, independently of the accretion rate and viscosity coefficient. 



Introduction. — Coalescing binaries of supermassive 
black holes (SMBHs) are the primary sources of gravita- 
tional waves (GWs) for the planned Laser Interferometric 
Space Antenna (LISA [38(). Recent advances in numeri- 
cal relativity enable to forecast the GW luminosity and 
waveform during the final phase of a SMBH coalescence 
event (e.g. Ref. [j], 0, Q and references therein). The 
peak GW luminosity , Lqw ~ 10 56-57 erg s _1 , will be 
observable with LISA out to high cosmological redshifts. 

As the GWs propagate away from their source, they in- 
teract with matter in several ways. The shear they induce 
in surrounding gas can be dissipated through viscosity 
The GWs could also drive transverse and longitu- 
dinal density waves H, excite resonant vibration modes, 
boost the frequency of photons 0, 0] , lead to graviton- 
photon conversion [8, 9], and couple to Alfven and mag- 
netohydrodynamic waves in strongly magnetized plasmas 
13, EH, E lH, El • These interactions are so weak that 
GWs are expected to escape from the densest environ- 
ments like the cores of supernova explosions, gamma 
ray bursts, or the early universe, and travel across cos- 
mological distances without any noticeable attenuation. 
However, in the vicinity of coalescing SMBH binaries, 
even a miniscule coupling with matter could lead to a 
bright electromagnetic (EM) signal. In this Letter, we 
demonstrate that the viscous dissipation of GWs in the 
astrophysical environments of SMBH binaries might be 
detectable. 

The merger dynamics of a pair of gas- rich galaxies with 
SMBHs generically channels large quantities of gas to the 
central region and creates a gaseous envelope around the 
resulting SMBH binary 15]. The presence of gas and 



stars is expected to catalyze the hardening of the SMBH 
binary 



17[ to a separation where GW emission 



on its own - would be capable of shrinking the orbit 



on a Hubble time [18(. The high infall rate of gas in 



this environment is expected to lead to the formation of 
a geometrically-thin accretion disk 17]. For binaries of 



nearly equal mass SMBHs, the tidal field of the binary 
would open a central cavity in the disk with a radius of 
about twice the orbital radius of the binary [19] . During 
the final phase of SMBH coalescence, the cavity would 



not be able to track the ever increasing rate at which 
GW emission shrinks the binary orbit and so the cavity 
radius would freeze at r m ; n ~ 12Qrsq$ A5 M® m . Here, 
rs = 2GM/c 2 , q = 4g/(l + q) 2 is the scaled symmetric 
mass ratio, and M7 = M/(10 7 M Q ) is the total binary 
mass. The GW dissipation in this punctured disk should 
inevitably lead to an EM transient. 

If the prompt EM counterpart to a GW signal is suf- 
ficiently bright, it would enable observations of SMBH 
mergers with traditional telescopes long before LISA be- 
comes operational. Its detection would provide a unique 
test of general relativity in the strong-field regime. A 
successful identification of an EM counterpart to a LISA 
source would have far-reaching consequences [2(i| . 

Recent studies considered other mechanisms for EM 
counterparts to SMBH mergers: (i) periodic variation in 
the gravitational potential owing to the orbital motion of 
the binary during the early stages of the inspiral 21] ; (ii) 
shocks induced by the sudden mass loss of the binary due 
the final GW burst 22}; (Hi) shocks induced by a super- 
sonic gravitational recoil kick I23l. \24 , \25 \; and (iv) infall 
of gas onto the SMBH remnant [19] . We expect the vis- 
cous GW heating effect to dominate over mechanism (i) 
during the late inspiral at distances much larger than the 
binary separation, since the tidal potential driving effect 
(i) scales as cP/r 3 . The viscous dissipation of GWs is 
unique in its ability to yield a prompt EM counterpart 
within hours-days after the peak GW burst, because it 
is driven by the time-evolution of the GW luminosity. 
In contrast, the mass loss effect of mechanism (ii) re- 
quires the orbital timescale of days-weeks necessary for 
the gas to respond and shocks to build up 24]. The 
supersonic kick from mechanism (Hi) becomes effective 
after the disturbance to the disk propagates out to the 
radius where it is supersonic, over a timescale of months 
years [lH, 0, IIH . Finally, infall of gas in mechanism (iv) 
occurs only years after coalescence (l9] . 

Gravitational waves from a black hole merger. — The 
GWs produced by a black hole binary merger are dom- 
inated by the spin-2 weighted I = 2, m = ±2 spherical 
tensor harmonic 0, 0] • We approximate the GW energy 
flux by its asymptotic behavior at large distances [26| . In 
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spherical coordinates, 



eow(t,r,M) = Y (0) 



rci 

4 7rcr 2 



(i) 



where 9 is the angle relative to the total angular mo- 
mentum vector, Y{0) = (5/2)[sin 8 (6»/2) + cos 8 (6>/2)] has 
a unit average on the sphere, and icw^rct) ^ s ^ ne GW 
luminosity at retarded time t re t = t — r/c, which we ap- 
proximate as Lgw(£) oc god^o* — ^ i I i ) " 5 / 4 for t < 0, 
and exp(— ci/(2.5rs)) for i > £i, while being constant in 
between. We set the overall scaling and t\ to match the 
Newtonian inspiral luminosity for ( « 0, and the peak 
luminosity from numerical simulations at t 
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spectively (see Fig.[2]below). This is correct within a fac- 
tor of ~ 2 depending on the magnitude and orientation 
of SMBH spins @, y| . The total GW energy released is 
AP GW = L Gw dt = kMc 2 , where 3% < n/ql < 7%, 
and 5% using our simple fit. The timescale of the intense 
GW burst is Ai GW = nMc 2 /L™ k ~ 20r 5 /c. 

GW dissipation in a viscous medium. — The stress- 
energy tensor of a viscous medium is augmented by 
T^ v = —2-na^, where rj is the dynamical shear viscos- 
ity coefficient and is the fluid's rate of shear [27|. As 
GWs traverse the fluid they induce a shear <x M y = \h^ v 
[39l ] . where an overdot denotes a time derivative of the 
metric perturbation, h^. The weak- field Einstein equa- 
tion yields Dh^ — — WirCr/h^fc 4 0, [53|, implying that 
the GW energy density ecw is dissipated at the rate, 
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eheat = — eGW = n — ecw- 



(2) 



Thus, the energy density dissipates exponentially with a 
time constant of td — (IQnGrj / c 2 )^ 1 . 

Heating a Thin Accretion Disk. — In the standard 
"a-model" of radiatively-efficient accretion flows 2^, 29|, 
|30| , the gas orbits around the central SMBH within a thin 
co-planar disk, characterized by a vertical height H(r) <C 
r and a low temperature T < 10 6 K. The viscosity coef- 
ficient is parameterized as rj(r) = (2/3)oP(r)/f2(r) [ioj . 
where 0,(r) = GM/r 3 is the angular velocity, a < 1, 
and P(r) is either the gas pressure, P gas , or the to- 
tal (gas+radiation) pressure, Ptot- We therefore write 
P(r) = Pg as P t 1 7 fc with < b < 1. The mass accretion 

rate M can be expressed in terms of the Eddington lu- 
minosity Le(M) as M = ttiMe, where Me = eLg/c 2 
and e is the fraction of the accreted mass which gets con- 
verted into radiation. Once the disk opacity k is speci- 
fied, the accretion model is fixed by the free parameters 
to, e, 6, and M 31, 3^]. The physical characteristics of 
the gas, such as the midplane temperature T(r), surface 
mass density S(r), and the scale-height H(r), can be all 
expressed in terms of these parameters. 

The GW energy absorbed per unit area of the disk sur- 
face, ifeheat, depends on the combination r/H in Eq. [[2"]). 



Remarkably, this particular combination is simply related 
to the mass accretion rate M in a steady state [3fJ , 



r)(r)H(r) 



i/(r)E(r) _ M 
2 ~ 6tt 



= const, 



(3) 



where, v = rj/p = 2Hr//T, is the kinematic viscosity and 
p is the mass density of the gas. 

Substituting Eqs. ^ and {3} into Eq. ©, we find that 
the rate of GW dissipation per unit surface area is inde- 
pendent of the disk viscosity or opacity, 
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(4) 



Thus, equal amounts of heat are dissipated by the GWs 
per logarithmical radius bin of the disk. The heating rate 
is proportional to M but is otherwise independent of the 
disk parameters. 

The heating rate should be compared to the stan- 
dard dissipation rate of the accretion disk, Pedisk = 
(3/87r)GMM/r 3 , 



^heat (^r 



edisk(V) 
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Y(d)r 3 L^(t ret ), 



(5) 



where 



£gw/(10- 3 c 5 /G). In general. 



t GW 



1 



during the peak emission of GWs, independently of M. 
Thus, the peak GW heating rate exceeds the standard 
disk dissipation rate at r$ > 1 by a factor of (32/9)!>3. 
If the disk resides in the orbital plane of the binary, 
ehcat/edisk = (10/9)r3< This result is universal and inde- 
pendent of M or the accretion disk parameters. 

Cooling Timescale. — The excess heat deposited by 
GWs will eventually be radiated away electromagneti- 
cally. The EM light curve depends on the uncertain de- 
tails of the turbulent accretion disk, and in particular on 
the vertical transport of heat. 

For an optically-thick disk with a vertical optical depth 
r(r) ^> 1, the diffusion timescale of photons out of the 
midplane is t^s ~ tH/c [111, [33| . The timescale for a 
patch of the disk to change its thermal energy content by 
turbulent heat transport is ftherm ~ idyn/o, where idyn ~ 
H/c s ~ l/Q and c s is the sound speed [34|. The excess 
dissipation of heat by GWs will eventually be radiated 
away within a timescale, t c = min(idiffj ^therm) 13] 41]. 

We calculate the diffusion timescale for various disk 
models, following Refs. 31, 3^]. We consider a viscosity 



term that is proportional to the gas or total pressure, 
an opacity coefficient that is dominated by electron scat- 
tering with k cs = 0.35cm 2 g _1 , or free-free transitions 
with Kg = 6.9 x 10 22 pT~ 7 / 2 [cgs], and Eddington accre- 
tion rates of to = 0.1 and 1, respectively. We also adopt 
a = 0.3 and Mj = 1 or 10. Figure Q] shows the result- 
ing tdis(r). The disk can be divided into three distinct 
regions: (a) an inner region dominated by radiation pres- 
sure; (b) a middle region dominated by gas pressure and 
k cs ; and (c) an outer region dominated by gas pressure 
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FIG. 1: The radiative diffusion time, idifi, over the geomet- 
ric light travel time, r/c, as a function of radius for various 
a-models with a = 0.3. The solid and dashed curves corre- 
spond to rh — 0.1 and 1, respectively. Two sets of curves are 
shown for Mr = 1 and 10, respectively, with the lower mass 
corresponding to the lower curves. 



and Kg. In region (a), the viscosity prescription affects 
the results dramatically. If r\ oc P gas then tdm ~ ^dyn, 
whereas if rj oc P to t then radiative diffusion is inefficient 
and t c ~ itherm- For typical parameters, the outer re- 
gion corresponds to very large radii, at which the disk 
may fragment due to its self-gravity. We find that for 
to = 0.1 and 0.1 < M 7 < 10, the range of 

3 < ct c (r)/r < 30, (6) 

applies to all radii of interest. We note that the cooling 
time might be reduced if the radiation escapes through 
regions of low gas density, if a significant fraction of the 
turbulent magnetic energy is dissipated in surface layers 
[3o| | , or if a photon-bubble instability operates [36| . 

Brightening of the Disk. — For stationary disks, there 
is no radial heat transport on the relevant timescales, 
since the viscous time is long, t c <C t v - lsc ~ r 2 jv ~ 
tthevm(r/H) 2 The flux excess due to GW heat- 

ing, AF(t,r) = F(t,r) — F<u s k(t,r), is determined by an 
ordinary first-order linear differential equation in time t, 

t c AF + AF = He heat , (7) 

where the driving term is determined by the instanta- 
neous GW luminosity, Lgw(*) ; from Eq. ([4]). The ob- 
served flux depends on the disk parameters through t c (r). 
The results are particularly simple in two limiting cases, 
namely when t c is much smaller or much larger than 
the heating timescale, for which AF(£ ret , r^/F^iskM = 
ehoat/edisk x {1 or At GW /t c exp(-i rot /i c )}, respectively. 



FIG. 2: The excess luminosity curve before (t < 0) and after 
(i > 0) the binary coalescence event, relative to the luminosity 
of the punctured disk. The time axis is shown on a logarithmic 
scale at both negative and positive values (in units of rs at 
the bottom or 'M? days' at the top). The solid and dashed 
curves correspond to cases where the accretion disk is face on 
or edge on, respectively. The top curve applies if the cooling 
time is shorter than the geometrical delay r/c at each radius, 
while the bottom curve show the case where disk cooling is 
limited by photon diffusion. The dotted curve depicts L°¥(t) 
for reference. The j/-axis is normalized by the luminosity of 
the punctured disk for a cavity radius, r m t n = 10 2 r m i nj 2?"s. 

Substituting Eq. (O in these simple cases gives 

AF(< rct ,r) ^ f f Y(6)r 3 L G f (t ret ) if Ai GW » t c 
.FdiskM ~ \ f K Y(6)r/(ct c ) if A£ GW < t c 

(8) 

For example, substituting = 1, Y = 1, k = 3%, 

and r/(ct c ) = 0.1, we find that AF(r)/F disk (r) - 3.6r 3 
and 0.005 in the two cases, respectively. Depending on 
the binary orientation relative to the disk, Y, the result 
can be a factor of ~ 3 larger or smaller (see Eq. [1} . The 
net excess in the apparent luminosity of the disk is 

AL(t) = cosJW J J Ai?(4t)r)rd0d . skdn (g) 

r min 

Here, t' Iet = t - (r/c)(l - sm6» obs cos <?!) d i sk ), where 9 ohs 
is the inclination of the disk relative to the line of sight, 
0disk is the azimuthal angle within the disk, and d^z) is 
the luminosity distance to the source at redshift z. 

We integrate Eq. J7]) exactly for our fit to LqwW an d 
find the luminosity excess AL{t) using Eq. ©. Fig- 
ure [2] shows the result relative to the punctured disk lu- 
minosity Ldisk- We assume that the disk has an inner 
cutoff at 10 2 r m i nj 2rs, and adopt (qo, to, a, b, M7, Y) = 
(1, 0.1, 0.3, 0, 1, 2.5), We consider two cases: (i) t c <C r/c; 
and (ii) t c {r) for the a-disk model. The dotted (blue) 
and solid (red) curves correspond to o bs = and 7r/2, 
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respectively. The light -curve is highly sensitive to # bs 
due to the geometric (GW travel-time) delay. If observed 
in the disk plane, the disk flare at all radii is seen co- 
incidently along the line-of-sight to the SMBHs, but 
with a delay in other directions. However, if the disk 
is observed face-on, then different annuli are seen coin- 
cidently and the geometric delay increases with radius. 
In the latter case, the peak EM luminosity is delayed by 
ic(^min) + (?*min/c) cos $ bs relative to the GW peak. The 
light curve has a characteristic shape, with AL oc |i| -5 / 4 
in the inspiral phase, AL oc i _1 after the peak for a time 
~ ^cooi ~ 10r max /c, and an exponential decay at later 
times. For unequal masses, AL is reduced by q$. 

Discussion. — Figure [2] implies that the excess lu- 
minosity of a thin circumbinary disk peaks with a de- 
lay ~ IOM7 hours relative to the peak in the GW 
burst. The magnitude of the excess luminosity peak 
is AL/Ldisk ~ (fmin/10 3 rs) for t c -C r/c, and a fac- 
tor of - 130to 4 / 5 M 7 1/5 smaller for more realistic t c (r) 
values. This amounts to ~ (10 _4 -10~ 3 ) x During 
the i _1 decline of the EM transient, the characteristic 
emission wavelength corresponding to the surface tem- 
perature of the disk ~ 3.6 x 10 3 r 3 3 ^ 4 K is in the infrared 



band A ~ l.Sr 3 / 4 /im, and increases in proportion to t 3 
as the GW propagates outwards. Future monitoring sur- 
veys, such as PAN-STARRS or LSST 0, could search 
for the expected transients [201 ] . 

In radiatively-inefficient (geometrically-thick) accre- 
tion flows [37[, the heating rate is also given by a uni- 
versal expression similar to Eq. ([5]). However, due to the 
low radiation efficiency, the resulting light curve is fainter 
and difficult to observe. We also find that GW heating 
is not sufficient for reversing the flow away from the BH. 

Measurement of the GW heating effect would provide 
an indirect detection of GWs with traditional electromag- 
netic observatories, and test the theory of general rela- 
tivity for the interaction of GWs with matter. Identifica- 
tion of an unambiguous EM counterpart to a LISA mea- 
surement of GWs, would enable to determine the source 
redshift and location as well as to constrain the gaseous 
environments of merging SMBH binaries. 
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